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| 

{ This report docunents the results of work accomplished by 
Systems Division, Vought Corporation, P. 0. Box 5907, Dallas, Texas 
75222, under Contract No. F08635-76-C-0237 with the Air Force Arma- 
ment Laboratory, Armament Development and Test Center, Eglin Air Force 
Base, Florida. Lt. Vayl S. Oxford (OLJC) managed the program for the 
Armament Laboratory. This effort was conducted during the period 

from May 1976 to October 1976, 
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: | publication. 
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[ SECTION I 
INTRODUCTION 


of 
This report and user's guide presents a computerized version 
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of the wing-mounted single store airloads prediction methodology found 
in External Store Airloads Prediction Technique (Reference 1). The 

ea prediction technique reported in Reference 1 is somewhat time-consuming 
and tedious to use in its present form due to the multitude of manual 


computations required to complete a full six-component airloads solu- 


tion. The computer routine developed during the current program and 
: reported herein will eliminate most of the engineering effort required 
| to evaluate captive store aerodynamic loads not only for current pro- 
duction aircraft-store configurations but also for new aircraft-weapons 
designs without resorting to expensive wind tunnel or flight tests. 
nS The technique should prove especially valuable for evaluation of new 

aircraft-weapons design in the preliminary design stages because of 

| the rapid response now available using the computer program for chord- 
wise, spanwise, and vertical variations in store location. Therefore, 


it permits trade studies to be conducted to minimize installed loads, 


for example, or to determine a location where the captive loads pro- 
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mote favorable separation. + 


The yawing moment method of the basic technique has been 
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extended during the current effort to include a broader range of angle 


of attack than previously presented. The improved yawing moment pre- 


diction has been incorporated into the computer program. 


The computer routine developed during this program is 
an accurate representation of the manual technique Plus improvements 
which includes $ix-component Single carriage captive store airloads 
Prediction capabilities for the basic airload as well as the incre- 


mental airloads due to aircraft yaw and adjacent store interference 


as detailed in Reference 1, 
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SECTION II 
COMPUTERIZATION OF STORE AIRLOADS PREDICTION TECHNIQUE 


A computer program has been developed that permits rapid 
prediction of store aerodynamic loads in the captive flow environ- 
ment for wing-mounted single stores requiring only a knowledge of 
store isolated aerodynamics, flight conditions, and aircraft-store 
geometry. 

The computer program was designed by considering the general 
and specific technical requirements outline by the Armament Jevelop- 
ment and Test Center (ADTC) and molding these requirements into a 
routine that is simple to input, duplicates the manual technique, and 
has the minimum possible storage requirements and execution time. 
GENERAL REQUIREMENTS 

General requirements for the subject effort were to develop 
a computer program that permits use of the wing-mounted single store 
portion delineated in AFATL-TR-75-87 (Reference 1) on the CDC 6600 com- 
puter at ADTC and, in addition, to extend the yawing moment method to 
include the full -4 to +12-degree angle-of-attack range included by 
the other five components. These requirements have been accomplished 
and are reflected in the computer routine described in this report. 
APPROACH TO COMPUTERI ZATION 

Because the basic technique involves many steps of compu- 
tation and data handling, this situation presents many options as 
to the formulation of a computer program. Some of the pertinent 


factors considered are: 


aerate tn 


1. Compromise between versatility and ease of input. 

2. Readability with manual technique (Reference 1). 

3.  Qutput options to meet military and industry needs. 

4, er Mtive diagnostics to aid the user, 

5. Minimize storage requirements and execution time. 

The contractor's experience in aircraft/store design and 
certification programs aided significantly in evaluating the various 
options and in formulating an effective and efficient computer program. 

The overall objectives were considered best satisfied by a 
computer program consisting of a main controlling routine that 
calls various subroutines to predict the airload components and 
provide supporting methodology. The details of the program, including 
a simplified flow chart schematic, nomenclature, input requirements, 
output description, and a sample problem, are presented and discussed 
in Appendix A; however, some of the salient features of the input 


deck, multiple run capabilities, and output features are discussed 


here, 


Design of the Input Stream 

The input data requirements for the computer program 
include Mach number, store reference length, dimensions of the 
store, dimensions of the aircraft, location of the store on the 
aircraft, and the isolated store aerodynamics, Some other parameters 
required for routine execution, such as store refererice area, are 
calculated internally. 

In order to provide the user with the most gefc72] and 


useful output, the predicted airloads are presented as coefficients 


in slope and intercept form, as discussed later. This approach 
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allows deletion or alrcrart angl@ or attack, angle of sideslip, and 
altitude from the input stream, and allows the user to investigate 
the airloads at any number of flight conditions after obtaining the 
computer output. 

The input data requirements were divided into six basic 
groups as shown in Figure 1. Each card in the input stream carries 
one of the descriptive names of Figure 1 in the first few columns 


of the caru. A brief description of each of the six input groups 


follows: 
TITLE 


TITLE cards contain alphanumeric descriptive information 
pertaining to the aircraft-store configuration to be run. 


Up to four TITLE cards may be included. 
STORE 


STORE cards contain most of the required store geometry 
and isolated aerodynamic information. Up to five STORE 


cards can be required for a configuration, STORE] through 
STORES. 
A/C 


A/C cards contain all parent aircraft geometry requirements, 
Two A/C cards are required for each configuration, A/C] 


and A/C2. 


The AERO card contains the Mach numbers at which airloads 
predictions are required. Up to seven Mach numbers may 
be specified between Mach 0.5 and 2.0, 
OPTION - The OPTION. card is used to specify the components for 
which predictions are required. A second field on this 
card is used to specify the components for which the 
internally looked-up variables are to be printed. 
SPA, PPA - SPA. and. PPA cards contain segmented store projected areas : 


used in the initial prediction summation procedure. 


TITLE 


STORE 


A/C 


AERO 


OPTION 


SPA, PPA 


Figure 1. 


- Header information 


- Store parameters 


- Aircraft geometry 


- Mach numbers 


- Computer routine control 


- Store projected areas 


Organization of the Input Data 
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Qutput Features 


The output features include four main sections: input data 
listing, store airloads predictions, diagnostic messages, and optional 
print of all variables looked up in the data base. These four cate- 
gories of output appear in printed form, as shown in Figure 2. 

The input data playback is always included in the output 


listing for user convenience. With the input data attached to the 


tabulated data, guesswork is eliminated as to what case the predicted 
loads data applies. In addition, errors in the input deck are 
readily traceable through the input deck listing. 

The captive store airloads predictions follow the input 
data playback. The basic airload. plus the incremental airjoads due 
to aircraft yaw and adjacent store interference are grouped. together 
for each of the six components, side force, yawing moment, normal 
force, pitching moment, axial: force, and rolling moment. The six 
components are organized and printed in the identical manner found 
in the table of contents of Volume II, Book 2 of Reference 1. This 
organization aids in the readability and rapid association of the 
predicted values with the various sections of the manual] technique. 
The coefficient data are presented in column form for each Mach 
number requested. After all component coefficients are presented, 
Gistance from store center-of-gravity to forward lug, store 
reference area, and store reference length data are presented. An 
asterisk to the immediate right of a coefficient indicates. that 
the prediction was made using-.a-parameter from the data base that 
was out of range or beyond the span of data used to derive the 


basic technique. An out-of-range warning is issued in the diagnostics 


section. 


INPUT DATA PLAYBACK 


AIRLOADS PREDICTIONS 


DIAGNOSTICS 


‘OPTIONAL OUTPUT 


Figure 2, Output Structure 
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The next section of the output centains diagnostic messages 
for the benefit of the user. The diagnostics are divided into two 
types: general and out-of-range warnings. The general diagnostics 
include messages primarily about input errors. These errors can be 
either fatal or non-fatal. A non-fatal input error might be that 
two AERO cards were input and program execution was performed using 
the last card read. A fatal error might be that the store diameter 
field was left blank and, therefore, assumed zero. In case of a 
fatal error, the input data are printed along with the fatal error 
message(s) and:-execution is terminated. The out-of-range warnings 
are designed: to- inform the user that a prediction has.:been made with 
a parameter that is beyond the bounds of the empirical -data base used 

to develop the basic technique. Thus, the prediction:has. been made 

with an extrapolated parameter. The warning issued:-contains the 

appropriate figure number from Reference 1 and prints the x, y value 

used in the prediction from that figure. It is then--up: to the user 

to consult the figure and apply his engineering judgement to decide 

if use of the extrapolated data overly compromises the final predictions. 
The: final -portion of the output stream contains the optional 

output printedi-only if requested on the OPTION card. The optional 

output contains..al] parameters looked up from the data: base for each 

airload. component computed. The optional output wil print the 


figure number and: x, y value used in the computation: from the data 


base (most figures in Volume II, Book 2, Reference 1). 


Multiple Run Capabilities 


Up-:to five separate aircraft~store configurations can be 


run in a single input deck by separating the stacked cases with END 
CASE cards as defiiied in Appendix A. If the same aircraft-store 
combination is involved, then only those input cards affected by the 
change need to be redefined. Thus, the input is constructed in a 


manner similar to Namelist input in FORTRAN, A maximum of seven Mach 


numbers can be run for each of the five cases. 


SECTION III 
APPLICATION OF COMPUTED RESULTS 


The tabulated coefficient results of the computer program 
present the basic airload (i.e., the captive store airload generated 
by a zero-yaw pitch excursion of the parent aircraft) and the incre- 
mental airloads due to aircraft yaw and adjacent store interference. 
The incremental airloads due to yaw and interference are predicted 
as increments to be added to the basic airload. Although the method 


of combining these effects was detailed in Reference 1, it is repeated 


here for convenience. Therefore, at a particular Mach-number the total 


captive airload experienced by..a store can be obtained from the 
foltowing generalized coefficient expression: 


Cc, #0, Fal, + Be + AC, (1) 


TOTAL —_ BASIC 8s INTE 
where: 
X - Can be y, n, N, M, A, or & representing side 
force, yawing moment, normal force, pitching 


moment. axial force, and roiling-moment, respec- 


tively 
Cy - Basic captive airload generated by a zero yaw 
BASIC pitch excursion of the parent aircraft 
Os - Incremental airload due to aircraft yaw per 
S degree store yaw angle, Bo. 
Be - Store yaw angle equal to YA/C for a right-wing 


store installation and “Yare for a left-wing 


store installation 


ac, - Incremental airload due to the effect of adjacent 
INTF store interference. 
, The generalized coefficient equation presented above, 


Equation (1:), applies to both the slope and. intercept predictions 


made by the computer program. 
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SECTION IV 
MODIFICATION OF THE YAWNING MOMENT TECHNIQUE 


A typical variation of side force and yawing moment for 
a wing-mounted single store is shown in Figure 3 as a function of 
angle of attack. This type of yawing moment variation extends from 
subsonic flow up through the transonic regime for the experimental 
data available during the technique development reported in Reference 1. 
The extreme change in slope of the yawing moment as a function of angle 
of attack occurs in the 4 to 6-degree angle-of-attack range for most 
configurations. During the technique development phase of the previous 
program, the yawing moment prediction technique was developed: for the 
~4 to +6-degree angle-of -attack range. One objective of the current 
program was to extend the capability of the yawing moment method to 
apply from -4 to +12 degrees angle-of-attack. This has been accomplished: 
by adding to the previous technique. The previous method has been: 
retained and still applies in the -4 to approximately 6-degree angle- 
of-attack range. 

Extension of the effective angle-of-attack range for 
the prediction was accomplished by predicting the angle of attack at 
which the break in yawing moment occurs, and then correlating. the 
yawing moment slope for the 6 to 12-degree range using: techniques. 
established in the previous study (Reference 1). The resulting: tech- 
nique for the entire -4 to +12-degree angle-of-attack range consists 
of a predicted slope and intercept for the lower angle-of-attack ‘range 


(-4 to 6 degrees) and a predicted slope and intercept for the higher 


range (6 to 12 degrees). The actual range of applicability for each 
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Figure 3. Typical Captive Store Side Force and Yawing Moment 
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of the equations is determined by the predicted angle of attack at 
which the yawing: moment slope break occurs. 
This approach was adopted after studying the available 

experimental data for high- and low-wing aircraft. It was observed 

for high-wing aircraft that the yawing moment slope break occurred 

at the angle of attack where side force was zero, as shown in Figure 3. 
This loading indicates that the store running load distribution is 
producing a couple at this angle of attack with maximum yawing moment 
occurring when side force is zero. This relationship -holds true for 
the available experimental data for high-wing aircraft for Mach numbers 
from 0.5 to 1.0 at which point the extreme non-linearity disappears 
from the yawing-moment component. For low-wing aircraft, this simple 
relationship did not hold true; nevertheless, a relationship was 
discovered between the angle of attack where side force is zero and the 
angle of attack where the break in yawing moment slope occurs. This 
relationship varied with Mach number over the range 0.5 to 1.2, above 
which the extreme non-linearity disappears for the yawing moment 
component for Jow-wing aircraft. The relationship between the angles 


of attack at-which the yawing moment slope break occurs -and the side 
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force is zero is presented in Figure 4 for both high- and: low-wing aircraft. 


With this relationship established for high- and low-wing 
aircraft, a method. was needed to apply the relationship- to other 
aircraft. This was accomplished by defining a wing-fuselage height 
parameter as shown in Figure 5. Armed with this definition, it is 
possible to determine for all wing-fuselage relationships the Mach 


number below which the yawing moment slope non-linearity occurs 


through the following equation: 
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Figure 5, Wing-Fuselage height Definition 
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If the Mach number for which an airloads computation is 
desired is less than MBREAK, the yawing moment non-linearity exists, 
and the procedure to ieee a new yawing moment slope and intercept 
that applies in the higher angle-of-attack range is required. This 
procedure is outlined in the following discussion. 

In order to determine the angle of attack where the yawing 
moment slope break occurs, it is necessary to know the angle of attack 
where side force is zero, Ps This parameter is determined from 


the following equation; 


(2) 
A ary: Oe al!) (3) 
Gee (=) 
q/a 
where; 
o*) - Basic side ferce intercept predicted from Reference 
q /a=0 
1 fer the Mach number of interest, ft?. 
(2) - Basic side force slope predicted from Reference 
Ge u 


ee : ft2 | 
1 for the Mach number of interest, deg 
Knowing ae from Equation (3), determine from Figure 4: 


a 
BREAK gy? deg 
WING 


a 
BREAK) oy » deg 
WING 


Then for the general case, 
h 


W 
a =o +— (a - ) + do (4) 
BREAK BREAK, oy Ne BREAK. roy BREAK, gy BREAK 
WING WING WING 
e 18 


Bw Eee 2 ee = os 
eee ote 
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22) Miscink hy eae 


aimee  sapeamitmdilitios| 
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where: 


hw - Wing-fuselage height parameter. Equal to 1.0 

ui for high-wing aircraft and 0 for low-wing 
configurations. Other configurations are defined 
by Figure 5, 


Adi oe aK - Final correction to BREAK to account for ae 


prediction tolerances (Figure 6), deg. 

The preceding equation for OppeAK determines the angle of 
attack at which the break in yawing moment slope occurs. For angles 
of attack less than Appear? the yawing moment siope and intercept 
prediction technique presented in Reference 1 is applicable. For 
angles of attack greater than Gppraxs the following equations predict 


the applicable slope and intercept: 


Yawing moment slope prediction for a > dporay: (see Figure 9) 


YM Nw 
my) stk (c K, ) * Kirpan] Sd + — ak (5) 
( ae SLOPE *°LOCAL Ay INTC he SLOPE tare 
where: 
Kot ope - Variation in an as a function of “Loca a,’ 
Figure 7; wae 
C LOCAL ~ Local wing chord of the parent aircraft at the 


location of the subject store, in. 


Ww : ce sin A 
hy Wing sweep correction factor, sin 45°? where A 


is the wing quarter-chord sweep angle, in degrees, 


of the parent aircraft. 
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Kite - Value of Ch where Cocarka = 0, Figure 7, 


a 7] eg 
S - Store saicbnks Aiea wd? ft? 
’ “7? 
d - Store diameter, ft 
tw - Wing-fuselage height relationship, previously 
Ne 


defined (Figure 5) 
AK - Incremental change in yawing moment slope due to 
the interference effect of the fuselage for high 


; ‘ : tt 
wing aircraft (Figure 8); deg. 


Yawing moment intercept for a > Cape AK: (see Figure 9) 


(Taro, : Co a, : (Fay! OBREAK * (Fax, (6) 


where: 


() - Variation of yawing moment with angle of attack 
a 
1 


for a ¢ Oppreax 2S predicted by Reference 1 


eae £t* 
basic airload, deg 


(@) - Variation of yawing moment with angle of attack 
a 
2 for a > AppeaK 4s predicted by Equation (5) 


3 
defined above, es 


OppeaK 7 Angle of attack where the yawing moment slope 
break occurs, Equation (4), deg 


Value of yawing moment at a=0 for a< OperaKk 


> 
-~ 
2 |= 
R 
n 
oO 
1 


as predicted by Reference 1 basic airload, ft°. 


The derivation of (“) 


0,° Equation (6), can be easily 

2 

understood by referring to Figure 9 from which the following expressions 
can be derived : 


20 


eee FT EE 
Steet eaeneenecemnie eae fi aaa eden enn meena 7 =: 


FIN SPA 


Curve 1: Applies if ame 


where FIN refers to the store WING for 
case ] stores and TAIL for case 2 stores 
of Reference 1. 


Curve 2: Applies if Lie > 82.0 in. 
Curve 3: Applies for all other configurations. 
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Figure 6, Variation of Alp pe ay with Mach Number 
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Figure 7; Yawing Moment Modification Parameters - K and K. 


SLOPE INTC 


22 


Pe? OD et AP ek Eee SLO EO MNS a Gee git Non Ot PE rs Py ee a. ak ae 
¢ 


| f=) 
ao 
| 1° 
| : 
f 1} = 
| ; = 
© & 
if : 3 — 
it ‘ i” 
ai 1 
b 1 3 c 
F ahi S 
[ } » 
: rs) 
| 3 
a 
x 
‘ ro) 
k oO 
- ct) 
© 4 Re 
‘s bal 
a) 
qe 
a 
; 7] 
f » 
é — 
ot 
t HH 
135 ov 
i < 
pi 
7 w 
! ” 
3 
rs 
i 
oS 
; @ 
: +S 
S 
Dm 
- 
| Le 
ft 
i 
i 
F 
: ‘ 
/ en tee 


9 Se ere ge er 


+e 


i madame 
en RE OR RE RT RT I TE 
os x 


© Experiment 


. ‘ Praere YM 
Figure 9. Guide to Derivation of (JD a=09 


24 


NT EE EL EF EE EEE TE WON 


4) = Cais crean* (a) 00, (7) 

a (=) ao “BREAK (8) 
Then, 

eee = dy + dy (9) 


Hence, the sum represented by Equation (9) yields the results 
presented in Equation (6) for (4) o=0.." 

The yawing moment prediction method presented in Reference 1] 
provides a predicted yawing moment slope and intercept for the angle-of- 
attack range -4° <a £ Oppray: Equations (5) and (6) presented here 
predict yawing moment slope and intercept, respectively, for the range 
pecan < o < 12°. Therefore, a combination of the two techniques yields 
a yawing moment prediction method that is valid for the -4° to +12° 
angle-of-attack range as required. This extension to the technique of 
Reference 1 has been incorporated in the computer program, and results 
from the program reflect the combination of these two techniques for 
the applicable Mach numbers. 

The figures required to predict the yawing moment slope and 
intercept for angles of attack greater than OBR EAK have been presented 
in this report as Figures 4, 6, 7, and 8. The figure numbers jin the 
manual technique (Reference 1) range from 1 to 909, In order to 
incorporate Figures 4, 6, 7, and 8 presented here into the digital 
data base for the computer routinz, they were renumbered beginning 


with Figure 910 as follows: 


EE ME 


DIGITAL DATA BASE FIGURE NUMBER 


Figure 4 910 (Low-Wing) 
911 (High-Wing) 


Figure 6 915 
Figure 7 912 (Kes ope) 

913 (Kiyre) 
Figure 8 914 (AKo) ope ) 


INTF 
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SECTION V 
EXRORS IN THE MANUAL TECHNIQUE 


Juring the commuter program development phase, several 
errors were discovered in the manual technique reported in Reference 
1. These errors were primarily plotting errors affecting four figures 
in Volume II, Rook 2, of Reference 1. These plotting errors involve 
Figures 26, 49, 53, and 59. The ordinate scale for Figure 53(a) 
and (b) should be negative instead of the positive scale shown. 
Figures 26, 49, and 59 are misplotted or poorly faired. The correct 
values for these three figures are presented in digital form, as used 
by the computer routine in Table 1. 

Besides the plotting errors, one typographical error 

PM 


affecting the equation for (FP ae0, on page 275 of Volume II, 


Book 2, of Reference 1 was found, The last factor appearing on the 
right side of this equation should read as follows: 

K PPA L 

Com n 


bed 


The definitions of KINTC, and KSLOPE, which follow on the same 
page should alse contain this factor. 
It is suggested that these corrections be noted in 


copies of Reference 1 to avoid further misuse. 
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Table 1. CORRECTIONS TO THE DATA BASE 
OF VOLUME II, BOOK 2 OF REFERENCE 1 


Figure 26 Figure 49 Figure 59 

x y x y X y 

0.0 0.0 0.0 0.0 0.0 0.16 

0.1 0.0 0.375 0.0 50.0 0.187 

0.2 0.0 0.390 0.52 100.0 0.227 

0.3 0.0 0.40 0.77 150.0 0.288 

0.4 -1.0x107° 0.4] 0.93 200.0 0.380 

0.5 -9.0x10"° 0.42 1.02 225.0 0.430 

0.6 -1.1x1079 0.45 1.20 250.0 0.445 

0.7 ~1.9x1079 0.475 1.29 275.0 0.448 

0.8 ~3,3x107° 0.50 153) 300.0 0.440 

0.9 -8.0x107° 0.65 1.40 350.0 0.414 

1.0 -1.25x1074 400.0 ‘0.378 

is) -1.6x1074 450.0 0.345 
Epos 1.2 -1.84x1074 500.0 0.325 
1.3. -2.04x1074 

1.4 -2.2x1074 

1.5 -2.31x1074 

ie -2.43x1074 


X is the independent variable 


Y is the dependent variable 
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SECTION VI 
CONCLUSIONS AND RECOMMENDATIONS 


A computer routine has been developed during this study 
that provides a rapid method to determine six-component captive 
store airloads for wing-mounted single carriage stores. The 
technique represented by the computer routine offers installed store 
airloads over a broad range of flight conditions which include Mach 
numbers from 0.5 to 2.0, angles of attack from -4 to +12 degrees, and 
aircraft yaw angles from -8 to +8 degrees. Airloads coefficients 
predicted by the routine are an accurate representation of the single 
carriage method reported in AFATL-TR-75-87 (Reference 1). As such, 


the computer routine allows inexpensive and rapid computation of 


installed store aerodynamic loads with accuracy sufficient for 


i 


preliminary design and evaluation purposes. Thus, the need for some 
wind tunnel and flight tests may be eliminated. 

Modifications of the yawing moment component to extend its 
usefulness above six degrees angle of attack were completed and 
incorporated into the computer routine. 


It is recommended that predictions from the method contained 


oo aed s x 
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in the computer routine be correlated with experimental data for 
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nies 


aircraft-store configurations not used in the method's derivation. 


soa ¢ 


rear a, ee: 


If significant differences are observed for specific configurations, 

: the technique should be updated to be applicable to new configurations. 
In addition, it is recommended that the remainder of the technique in 
Reference 1, captive airloads for multiple carriage stores, be 


developed into a computerized version for Air Force and industry use. 
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OPERATIONAL GUIDE 
Following is a list of conments which may be helpful to the 


user until he becomes familiar with the operation of this routine: 


1. Multiple cases can be run by separating the cases with "END 
CASE" cards, 

2. A maximum of five cases can be run at one time, 

3. A maximum of seven Mach numbers can be run for each case. 

4, Pata cards that are unchanged from case N to case N+? need 
not be included in the case N+] input deck, 

5, Four title cards should be used. 

6. Input fields for data that are not required for a particular 
case snould be left blank. (For example, if there are no 


interfering stores, the input fields describing interfer- 


“ ing store parameters should be left biank.) 
7. Title cards must be ordered as desired in the output, 
8. Cards other than title cards can appear in any order desired. 
9, If "x" or "+" is not input in column 10 of STORE3 card, 
then the pianform area of the wing/tail surface as -defined 
; for computation of Reference 2 must be input, as well as. 
; the sum of all wing and tail planform areas. 
; 10. For type 1 stores, "wing" refers to the lifting surfaces. 


lifting surfaces and "tail" to the aft set. 


i 
{ For type 2 stores, "wing" refers to the forward set of 
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11. For aircraft with fuselage side inlets, the ninimum distance 
from the store to the inlet as measured ti the aircraft 


yz plane must be input. For aircraft with centerline in- 


lets, a value of "400." should be input in columns 31-40 


of the A/C1 card. 


cee ry 


12. Each store segment containing SPA or PPA can have more than 


oe 


One area type in it. (For example, IN, 2N, 2B, 3B, 4B, 


bas "= 
derma te in a n men t Patt ee tl dancer a 


4W may appear to describe a particular store where N, 


waa ote 


B, W represent nose, body, and wing SPA or PPA, respectively.) 
13. The captive airloads for any components listed in columns 
41-60 of the OPTION card will be calcuiated automatically. 

An asterisk appears immediately following any predicted 

slope or intercept in the output listing which was cal- 
culated using extrapolated or out-of-range data from the 

data base (graphs of Reference 1). 

15, The maximum number of SPA or PPA cards that can be input is 


40. 
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i 16. The maximum number of store segments allowed is 35, 
: 17. If the Side Projected Area is not equal to the Plan Pro- 
jected Area for each segment type: 
a. For every SPA card there should be a corresponding 
PPA card, even if the Plan Projected Area for that segment 
type is zero. 
b. For every PPA card there should be a corresponding 
SPA card, even if the Side Projected Area for that 


segment type is zero. 


{ 
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18. For unfinned stores, a value of "90" must be input in 
columns 41-56 of the STORE2 card. 

19, Mach numbers may be input in any order on the AERO card, 

20. If the optional data base look-up printout for Figures 4 and 


"5 of Reference 1 are desired, "IP" must be input within 


columns 41-60 of the OPTION card. 
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FORTRAN NOMENCLATURE 
In an effort to provide the user the most convenient definition of 
FORTRAN symbols, the definitions of these symbols appear as comments in the 
listing of the computer routine. (See "Routine Listing" portion of this 
appendix.) By presenting the definitions in this manner, the location of 
each variable is identified without the preparation of a special section 


within this manual, thereby maximizing its simplicity. 
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INPUT INSTRUCTIONS 
This section provides a detailed description of the parameters required 


as input to the routine. 


eek rematch I til meester ett e 


CARD PARAMETER COLUMN/FORMAT DESCRIPTION 
(FORTRAN SYMBOL 
TITLE 11-80/7A10 General information describing 


the aircraft-store configuration 
to be investigated. There 
should be four "TITLE" cards. 


STORE] STORE TYPE 10/11 Identifies the type of store 
(ICASE) (Case 1, 2, or 3) to be studied. 
Case 1 ~ store with lifting sur- 
faces at aft end. Case 2 ~ 
store with both forward and aft 
lifting surfaces, Case 3 ~ un- 
finned store. 


y i ue " 
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wees 


Ky 11-20/F10.0 Ratio of lift of the store nose 
to lift of wing alone; required 

(AKN) for Case 1 stores, (See line 

80, Table 1, page 38 of Refer-~ 

. ence 2.) 

os ky 21-30/F10.0 Ratio of the body lift in the 

(W) presence of the wing to lift of 
(AKBW) the wing alone; required for 


Case 1 stores, (See line 48, 
Table 1, page 38 of Reference 2.) 


31-40/F10.0 Ratio of lift of the wing in the 
(B) presence of the body to lift of 
(AKWB ) the wing alone; required for 
Case 1 stores. (See line 47, 
Table 1, page 38 of Reference 2.) 


KN/B 41-50/F10.0 Required input for Case 2 stores. ‘ 
(AKNB) (See page 29 of Reference 1 for ; 
description, 

Kr7B 51-60/F10.0 | Required input for Case 2 stores. 
(AKTB) (See page 29 of Reference 1 for 


description. ) 
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(FORTRAN SYMBOL ) a 


Required input for Case 2 stores. 
(See page 30 of Reference 1 for 
description, ) 


Required input for Case 2 stores. 
(See page 30 of Reference 1 for 
description.) 


CARD PARAMETER COLUNN/FORMAT DESCRIPTION 
STORE] Ky ig 61-70/F10.0 
| (AKWBT) 
KINTF 71-80/F10.0 
(AKINTF) 
STORE2 Store Length,L  11-20/F10.0 


(STOREL) 


Store Diameter,d 21-30/F10.0 
(STORED) 


Nose Length,L 


Fane aL) A 31-40/F10.0 
NOS 


Wing/Tail Leading 41-50/F10.0 
Edge Sweep Angle, 


Nein 
(ALEFIN) 


Isolated Store 51-60/F10.0 
Lift Curve Slope, 


C 

Ly 

(CLAISO) 

Freestream Zero- 61-70/F10.0 
Lift Drag Coef- 

ficient,Cp 

(CDOISO) 


Store C.G. 71-80/F10.0 
(XCG) 


Length of the subject store, in. 


Diameter of the subject store, 
in. 


Nose length of the subject store, 
in., as defined by line 41, 
Table 1, page 38 of Reference 2. 


Sweep angle of the wing leading 
edge (Case 1 store) or Tail 
leading edge (Case 2 store) of 
the subject store, degrees. (A 
value of "90." should be input 
for unfinned stores.) 


Lift curve slope, C. , of the 
a 


subject store in freestream 
flow. Obtained from experimen- 
al data or by using the method 
presented in Reference 2, 


The Cy of the subject store in 
0 

freestream flow. Obtained from 

experimental data or by using 

the method presented in Refer- 

ence 3, 


The distance from the subject 
store nose to its center of 
gravity, in. 


Se eg ah et eee we ee, oe eee 


CARD PARAMETER 


tener ld cnet Daaeerime erties 


STORE3 +,x,1 or 2 
( ICONFG) 


ace nroan aan 


SPA=PPA or 
Blank 
(ISYM) 


| Wing Planform 
| Area 

i (WINGPA) 
4 

i 


Tail Planform 
Area 
(TAILPA) 


10/A) 


14-20/A7 


21-30/F10.0 


31-40/F10.0 


COLUMN/FORMAT DESCRIPTION 
(FORTRAN SYMBOL) 


Control parameter describing 
symmetry of store wing/tail. 

"4" . used if wing/tail are in 
the plus orientation when sub- 
ject store is in the carriage 
position, "x" ~ used if wing/ 
tail are in this orientation 
when the subject store is in the 
carriage position. 


"7" . used if wing/tail are 
"symmetric" as defined on page 
365 of Reference 1, ut not in 
"+" or "x" orientation. 


"o" . used if wing/tail are 
“unsymmetric" as defined on 
page 365 of Reference 1. 


If the side projected area of 
each segment is equal to the 
plan projected area use 
"SPA=PPA", In this case only 
SPA cards (as described later) 
must be input. If the side 
projected area of each segment 
is not equal to the plan pro- 
jected area, these columns 
should be left blank. Then 
both SPA and PPA cards must be 
input. There must be an equal 
number of SPA and PPA cards. 


If "1" or "2" is input in col- 
umn 10 of STORE3 card, the 
planform area of the wing as 
used in Reference 2 must be in- 
put here, ine, 


For Case 2 stores, if "1" or 
"2" is input in column 10 of 
the STORE3 card, the planform 
area of the tail as used in Re- 
perience 2 must be input here, 
inc, 


ae 
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CARD PARAMETER 


(SUMPA) 


STORE4 Longitudinal 
Overlap Distance 
of Outboard In- 
terfering Store, 


XINTF 
(XINTFO) 


Lateral Separa- 
tion Distance to 
Outboard Inter- 
fering Store, 


YONTF 
(YINTFO) 


Diameter of Out- 
board Interfer- 
ing Store, O NTE 


(DINTFO) 


Longitudinal 
Overlap Distance 
of Inboard Inter- 
fering Store, 


XINTE 
(XINTFI) 


Lateral Separa- 
tion Distance 
to Inboard In- 
terfering Store, 


YSNTF 
(YINTFI) 


Diameter of In- 
-poard Interfer- 
ing Store, OTE 


(DINTFI) 


COLUMN/FORMAT DESCRIPTION 
(FORTRAN SYMBOL ) 


STORE3 Sum Planform Area 41-50/F10.0 


11-20/F10.0 


21-30/F10.0 


31-40/F10.0 


41-50/F10.0 


51-60/F10.0 


61-70/F10.0 


If "1" or "2" is input in col- 
umn 10 of STORE3 card, the sum 
of all wing and tail planform 
areas must be input, iné. 


Distance from the nose of the 
subject store to the nose of 
the interfering store as mea- 
sured along the longitudinal 
axis of the subject store, in. 
(Positive if the subject store 
nose is forward of the inter- 
fering store nose. ) Required 
if an interfering store is out- 
board of the subject store. 


Minimum lateral clearance be- 
tween the subject store and the 
interfering store as measu’ed 
in the planview, in. Required 
if an interfering store is out- 
board of the subject store. 


Diameter of the interfering 
store, in. Required if an in- 
fering store is outboard of the 
subject store. 


Defined the same as “XINTFO". 
Required if an interfering 
store is inboard of the subject 
store. 


Defined the same as "YINTFO". 
Required if an interfering 
store is inboard of the subject 
store, 


Defined the same as 

"DINTFO", Required if an in- 
terfering store is inboard of 
the subject store. 


CARD 


| 


AERO 


A/C) 


PARAMETER 
FORTRAN SYMBOL 


STORES Length of Out- 


board Interfer- 


ing Store, LINTF 


(ALINFO) 


Length of Inboard 21-30/F10.0 
Interfering Store, 


LINTE 
(ALINEL) 


Lug Spacing 
(ZPLNSP) 


Store Midlug, 


Xue 
(SML) 


Segment Length 
(SEGL) 


y 
(Y1PRM) 


Mach Numbers, M 
(REALM(*) ) 


Local Chord 
Length, COCAL 


(CLOCAL) 
Pylon Midlug, 
x 


CML 
(PXDCML ) 
Wing Semi-span, 


(WSS) 


COLUMN/FORMAT DESCRIPTION 


11-20/F10.0 Length of the interfering store, 
in. Required if an interfering 
store is outboard of the subject 


store. 


Length of the interfering store, 
in. Required if an interfering 
store is inboard of the subject 
store. 


31-40/F10.0 Distance between suspension lugs 


on subject store, in. 
41-50/F10.0 Distance from the subject store 
nose to the store mid-lug lo- 
cation as measured along the 
store longitudinal axis, in. 
51-60/F10.0 Constant length segments of 
the subject store used to de- 
fine SPA and PPA, in. (user 
defined), 
61-70/F10.0 The minimum distance from the 
side of the fuselage to the 
Side of the subject store as 
measured in the planview, in. 
(see Figure A~1), 
11-80/7F10,0 Mach numbers at which captive 
store airloads are desired 
(0.5 <M < 2.0), 
11-20/F10.0 The aircraft local wing chord 
at the buttline corresponding 
to the subject store location, 
‘in. 
21-30/F10.0 Mid-lug location (fraction of 
aircraft local wing chord) of 
the pylon on which the subject 
store is carried. 


31-40/F10.0 Aircraft wing semi-span, in. 


(see Figure A-1). 


CARD 


A/t) 


A/C2 


SPA 


PARAMETER 
(FORTRAN SYMBOL 


Y 
(YBIG) 


y Buttline, Yep 
(YBL) 


(v1) 


Lambda, A 
(ALAMDA) 


Pyton Height, Z 
(ZPH) 


Store to Inlet 


Distance, LINLET 


(ALINLT(*) ) 
h,/be 


(HWDHF ) 


Segment Number 
(NSEG) 


COLUMN/FORMAT DESCRIPTION 


41-50/F10.9 


51-60/F10.0 


61-70/F10.0 
11-20/F10.0 
21-30/F10.0 


31-40/F10.0 


41-50/F10.0 


14-19/16 
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Distance from the side of the 
fuselage to the wing tip mea- 
sured in the planview, in. (see 
Figure A-1). 


Distance from aircraft center- 
line to pylon centerline (i.e., 
buttline of the pylon), in. 
(see Figure A-1). 


Distance from the side of the 
fuselage to the pylon center- 
line, in, (see Figure A-1). 


Quarter chord sweep angle of 
the parent aircraft wing, deg. 


Distance from the lower surface 
of the wing to the bottom of 
the pylon at the mid-lug point, 
in, 


Minimum distance from the side 
of the store to the fuselage 
mounted side inlets of the air- 
craft in the Y-Z plane, in. 
(see Figure A-1). (A value of 
"400." should be input for air- 
craft with centerline inlet.) 


hy = distance from the bottom 
of the fuselage to the wing 
plane, in., at the chord posi- 
tion of the pylon mid-lug. h¢ 
= height of the fuselage, in., 
at the chord position of the 
pylon mid-lug. hw/hf = 1.0. 
for a high-wing aircraft and 
hw/h¢ = 0.0 for a lowwing 
aircraft. 


Beginning at the store nose, 
the number associated with the 
constant length SPA segment. 


Ee oy SY oto bela ar aad 
3 ay 


Neh en 


or ty *, si 
oe 


CARD 


SPA 


PPA 


PARAMETER 
FORTRAN SYMBOL 


Segment Type 


20/A1 
(NSTYPE) 


Area 


21-30/F10.0 
(SPAT ) 


Segment Number 14-19/16 
(NSEG) 
Segment Type 20/Al 
(NSTYPE) 


Area 


21-30/F10.0 
(SPAT ) 


COLUMN/FORMAT DESCRIPTION 


Describes the 
each segment. 


type of area in 


N ~ nose area segment 
B ~ body area segment 
Wo. Wing area segment 
T ~ tail area segment 


Side projected area of segment 
NSEG, in¢, 


Same as for SPA. 


Same as for SPA, 


Plan Brogcctes area of segment 
NSEG, ins 


OPTION. SF,NF,YM,PM,AF,RM 11-27/7(A2,1x) @mponents appearing in these 


END CASE 


(INPTOP (*) ) 


IP, YM,PM,AF,RM, 
SF NF 
( INPGRA(*) ) 


41-60/7(A2 ,k) 


columns will be calculated 
Choose from 1-6 components and 
list in any order, 


The components appearing in 
these columns are those for 
which table look-up data are 
printed in the output. Any or 
all of these components may be 
listed in any order. 

IP = Initial Prediction. 


This card separates up to five 
cases for a single input deck. 


! HIGH WING WITH CENTERLINE INLET 


LOW WING WITH SIDE INLET 


Figure A-1. Input Description 
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OUTPUT DESCRIPTION 


The tabulated output consists of four distinct sections of 
information. First, a playback of the input data appears followed 
by the second section containing the predicted captive airloads in 
coefficient form. Following the predictions, a diagnostics section 
is included that contains informative statements concerning extra- 
polated data used in the predictions, as well as descriptions of 
non-fatal input errors. The final output section contains the values 
of parameters obtained from the data base used in the computation. 
This output section is optional. The output is structured in this 
manner for all cases except when a fatal input error is encountered. 
In case of a fatal error, the input data are printed along with the 
fatal error message(s) and execution is terminated. 

Figure A-2 presents the typical output arrangement of a success- 
ful run, showing the input playback, predicted captive airloads, 
diagnostic messages, and data base values. Each of the four output 
sections is discussed in greater detail in the following paragraphs. 

A. Input Data Playback 
A playback of the input data is furnished to provide 
the user with a convenient check of the input data 
and for future reference in associating predictions 
with the corresponding input request. 

8. Captive Airloads Predictions 


The final predicted coefficients. for the subject 


configuration are presented in column form for the 


LISTING OF INPUT DATA 


SUMMARY OF PREDICTED SLOPES 
AND INTERCEPTS FOR THE 
OUTPUT OPTION CHOSEN 


DIAGNOSTICS: 
(1) INPUT ERRORS 


(2) LIST OF OUT~OF-RANGE 
DATA USED 


— 


DISPLAY OF ALL VALUES 
OBTAINED FROM DATA. BASE 
USED TO CALCULATE THE 
CASE (OPTIONAL) 


Figure A-2. Schematic of Output 
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Mach numbers requested. The sign convention for 

the predicted airloads is identical to the single 

carriage sign convention presented in Volume II, 

Book 1 of Reference 1. Moments are referenced to , 
the store center of gravity. The coefficients are : 
presented as a basic slope and intercept for each 

component plus incremental slopes and intercepts 

due to the effects of aircraft yaw and adjacent 


store interference. 


Whenever extrapolated data from the data base are used 
in the calculation of a particular coefficient, an 
asterisk will appear immediately to the right of that 


coefficient. An exception to this procedure involves 


out-of-range data from Figures 4 or 5 of Reference 1. 


Ca 
. 


i 


Since these figures are used in the initial predictions 
portion of the computer routine, the flags are not 


issued. Predictions are flagged only if the out-of- 


ee 
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range parameter involves an empirically derived 


. correction from the experimental data base. 

¥ C. Diagnostic Messages 

i Various error messages have been incorporated into 
ql the routine to inform the user of any errors, both 

4 fatal and non-fatal, that exist in the input stream. 


attention to less serious mistakes in the input deck 
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| In addition, warning messages are included to call 


and to inform the user that a prediction has been 
made with a parameter that is beyond the bounds of 
the empirical data base used to develop the basic 
techinque. The out-of-range message applies to 
those predicted ceefficients that are flagged with 
an asterisk. 

Optional Output 

Th: presentation of the data base parameters used 
in calculation of the captive airloads is the last 
portion of the output. In addition to the "X" and © 
"Y" values from the data base graph, the appropriate 
figure numbers and line numbers are also presented. 
This portion of the output is optional and can be 
supressed by leaving columns 41-60 blank on the 


option card in the input deck. 


The figure numbers appearing in this portion of the 
output correspond to the figure numbers in Volume II, 
Book 2 of Reference 1, which forms the digital data 


base used by the routine. 


The line numbers presented here can refer to a variety 
of situations, depending on the type of graph from 
which data are obtained. Obviously, if a given 

graph in Reference 1 has more than one line, some 
distinction of data lines must be made. If there 

is only one line on a given graph, it is always 


labeled line 1. 


48 


ah Se a cr ie Se ree ee eS es 


— a ET 
rid eee eR EE Ow = 
Spe ye 


If the graph consists of several lines, each designating 
data for different Mach numbers, the lines are dis- 
tinguished by assigning line numbers in order of 
ascending Mach numbers. That is, line 7 will corres- 
pond to the lowest Mach number curve present, and 

the largest line number will correspond to the highest 
Mach number curve. For example, the data base for 


Figure 38 in Reference 1 will contain Lines 1 through 4. 


For graphs that contain incremental data due to air- 
craft yaw, the curve for negative store siaeslip (-8s) 
is line 1, and the curve for positive store sideslip 


(+8s) is line 2 as in Figure 95 of Reference 1. 


For graphs containing data describing the incremental 
NN load due to adjacent stove interference, the curve 
presenting data for the effect of the outboard inter- 
fering store is }ine 1. Line 2 is comprised of the 
data concerned with the effect of the inboard inrter- 
fering store. Figure 55 in Reference 1 illustrates 
typical curves defining the effect of outboard and 


: inboard interfering stores. 


If, as in Figure 48 of Reference 1, there is a graph 
with both gs and Mach number data, the line numbers 
increase with ascending Mach number, but there are 


two sets, one for (+8s) and one for (-gs). The set 


for (-8s) will always appear before the set for (+8s) 


SESS ee ee ae oe ee ee A 


a ES 


in the optional output. So, even though the output 


shows, for example, two "Figure 48, line 3" data sets, 


by the definition presented here the first corresponds 


to the (-gs) calculation end the second to the (+3s). 


The same type of situation can result when both inboard/ 

outboard interference and Mach number data appear on : 
the same graph, as in Figure 56 of Reference 1. In 

this case, the first set of data printed for a given 

figure number is for the calculation using the outboard 

interfering store curve and the second set is for the 

inboard. Obviously, both sets will appear only when 

there is an inboard and outboard interfering store 


for the case being computed. 


ROUTINE DESCRIPTION 


This section contains a brief description of the function of each 


portion of the computer routine. 


1. BLKRED 
BLOCK DATA routine stores values of parameters that do not 
change during the execution of the program. 

2. BLKFIL 
BLOCK DATA routine loads the file, identifying the variables 
with a file number. 

3. STORES 
External store airloads prediction technique main routine. 
STORES is the controlling program which cycles on up to five 


cases of data. 


SUBROUTINE DESCRIPTION 

1. CHECKR 
CHECKR assigns an input code to each card in the input stream. 
It checks for certain variables in which a blank field is valid 
and also checks for errors in card identification. 

2. READIN 
This routine reads data according to the format chosen by the 
input code that is assigned in CHECKR. It also checks fields 


for valid data and sorts SPA and PPA cards into ascending 


numerical order. 


ce 


neg 


INITAL 
INITAL sets up most constants used in the program as well as 
calculates the initial predictions for side force, normal 
force, pitching moment, and yawing moment. 

PRINTR 
PRINTR copies the files to output. 

ERROR 
This routine prints to the error file the appropriate message 
for each input error or warning encountered. 

LININT 
LININT is a linear interpolation scheme that retrieves values 
in the data base necessary to perform the required airloads 
computations. In addition, it contains the optional graph 
printout scheme. 

TWOD 
TWOD linearly interpolates between two graph lines in the data 
base for a specific independent variable. 

STOPRT 
This routine prints al] predictions for the specified Mach 
numbers containing the calculated output. 

SIDFRC 
SIDFRC calculates the predicted side force basic slope and 
intercept as well as the increments due to aircraft yaw and ‘ 
adjacent store interference for Mach numbers between 0.5 and 


2.0. 


10. YAWMOM 
YAWMOM calculates the predicted yawing moment basic slope and 
intercept, as well as the increments due to aircraft yaw and 
adjacent store interference for Mach numbers between 0.5 and 
2.0. 

11. NORFRC 
NORFRC routine calculates the predicted normal force basic 
slope and intercept, as well as the increments due to aircraft 
yaw for Mach numbers between 0.5 and 2.0. In addition, it 
calculates the incremental normal force slope and intercept 
due to adjacent store interference for a Mach number of 0.5. 

12. NORSBI 
This routine calculates the incremental normal force slope and 
intercept due to adjacent store interference for Mach numbers 
other than 0.5. 

13. PITMOM 
PITMOM calculates the predicted pitching moment basic slope and 
intercept as well as the increments due to aircraft yaw and 
adjacent store interference for Mach numbers between 0.5 and 
2.0. 

14, AXTFRC 

. This. routine calculates the predicted axial force: basic slope 


and intercept and the increments due to aircraft yaw and adjacent 


store interference for Mach numbers between 0.5 and 2.0. 


15. ROLMOM 
ROLMOM calculates the predjcted rolling moment basic slope 
and intercept and the incremental intercept due to aircraft 
yaw and adjacent store interference for Mach numbers from 0.5 
to 2.0. 

16. ALNLOC 
This routine looks up a dependent value from the data base 
assuming the independent values are evenly spaced. 

17. BREAKP 
BREAKP finds the Mach number break points above and below 


the subject Mach number for use in the applicable Mach 


number correction calculation. 
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INITIAL DISTRIBUTION 


Hq USAF/RDQRM 2 
Hq USAF/SAMI 1 
Hq USAF/XOXFM 1 
AFIS/INTA 1 
AFSC/DLCA 1 
AFSC/IGFG 1 
AFSC/SDZA 1 
AFAL/DHO 1 
AFWAL/Tech Lib/FL2802 1 
ASD/ENFEA 1 
FTD/PDXA-2 1 
AFWL/NSC 1 
AFWL/NSE 1 
AFWL/SUL 1 
AUL (AUL/LSE-70-239) 2 
DDC 2 
Ogden ALC/MMWM 2 
TAC/DRA 1 
57 FWW/DOS 1 
USAFTFWC/TA 1 
6510 ABG/SSD 1 
HQUSAFE/DOQ 1 
HQPACAF/DOO 1 
Naval Rsch Lab/Code 2627 1 
Naval Air Sys Comd/Code 530C 1 
Naval Air Sys Comd/Tech Lib 1 
Naval Surface Wpn Ctr/Tech Lib 2 
Naval Ordnance Stn/Tech Lib 1 
Naval Air Test Ctr/Tech Pubs 2 
USNWC (Code 533)/Tech Lib 1 
Sandia Lab/Tech Lib Div 3141 1 
Rand Corp/Library-D 1 
TACTEC 1 
USAFTAWC/TEFA 1 
TAWC/TRADOCLO 1 
ADTC/SES 1 
ADTC/SD23 1 
AFATL/DL 1 
AFATL/DLOSL 2 
AFATL/DLY 1 
AFATL/DLJ 1 
AFATL/DLJC 1 
AFATL/DLJF 1 
AFATL/DLJK 1 
AFATL/DLJM 1 
AFSC/SDTA 1 
ASD/ENESS i 
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